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Synopsis 

We introduce a model constitutive law for the dissipative stress tensor of lamellar 
phases to account for low frequency and long wavelength flows. Given the uniaxial 
symmetry of these phases, we argue that the stress tensor must be the same as that 
of a nematic but with the local order parameter being the slowly varying lamellar 
wavevector. This assumption leads to a dependence of the effective dynamic viscosity 
on orientation of the lamellar phase. We then consider a model configuration compris- 
ing a domain boundary separating laterally unbounded domains of so called parallel 
and perpendicularly oriented lamellae in a uniform, oscillatory, shear flow, and show 
that the configuration can be hydrodynamically unstable for the constitutive law cho- 
sen. It is argued that this instability and the secondary flows it creates can be used 
to infer a possible mechanism for orientation selection in shear experiments. 



I. INTRODUCTION 



Recent interest in block copolymers arises from their ability to self-assemble at the 

nanoscale through microphase separation and ordering, leading to mesophases with 

variou s types of symm etries, such as lamellar, cylindrical, or spherical Fredrickson and Bates! 
(1996); Larson ( 1999j )]. However, when processed by thermal quench or solvent cast- 
ing from an isotropic, disordered state, a macroscopic sample manifests itself as a 
polycrystalline configuration consisting of locally ordered but randomly oriented do- 
mains (or grains), with the presence of large amount of topological defects and un- 
usual rheological properties. The development of the equilibrium state characterized 
by macroscopic orientational order, as desired in most of applications, requires unre- 
alistically long times; hence external forces, such as steady or oscillatory shears are 
usually applied to accelerate domain coarsening and induce long range order. How- 
ever, the mechanisms responsible for the response of the copolymer microstructure 
and the selection of a particular orientation over a macroscopic scale are still poorly 
understood. In this paper we focus on the case of imposed oscillatory shear flows on 
lamellar phases of block copolymers, and present an orientation selection mechanism 
originating from an effective viscosity contrast between lamellar phases of different 
orientation. Our study is based on a mesoscopic or coarse-grained description of a 
copolymer, but our results are expected to apply to other systems with the same 
symmetry as a microphase separated block copolymer. 

The response of lamellar block copolymers to external shear flow can be classified 
according to three possible uniaxial orientations (see Fig. |TJ): parallel (with lamellar 
planes parallel to the shearing surface), perpendicular (with lamellae normal along the 
vorticity direction of the shear flow), and transverse (with lamellae normal directed 
along the shear) . Solid-like or elastic response is expected for lamellar phases of trans- 
verse orientation, while fluid-like or viscous response follows for the other two, leading 
to different rheologi c al properties at low shea r frequ encies as measured experimentally 
|Koppi et all (|l992f ); iFredrickson and Bated ( 1996f )]. It is known that both parallel 
and perpendicular alignments are favored over the transverse one under shear, and 
either of them would be ultimately selected by shear flow, as observed in most shear 
aligning experiments |Larson (|l999h ] (although some of e xperimental work indic ates a 
coexistence between parallel and transverse orientations Pinheiro et all (|l996[ )]. fact 
that might be a result of strong segregation and/or molecular entanglement). Of 
particular interest, and the least understood, is the selection between parallel and 
perpendic ular orientati o ns an d the dependence on shear fre quency u and strain am- 
plitude t iKoppi et al. (Il992l); Marine and Wiesnerl (1997)] as well as temperature 
|Koppi et all (|1992n ; IPinheiro and Winevl (|1998T )]. Near the order-disorder transition 
temperature Tqdt and at low shear frequencies, parallel alignment ha s been found 



inpoly(ethylene-propylene)-poly-(ethylethylene) (PEP-PEE) samples Ko ppi et al. 
( 1992f )]. while for pory(styrene) -pory-(isoprene) (PS-PI) co p olymers, th e obse rved 
ultimate orien tation is parallel Marine and Wiesneil (|l997l ): iLeist et ~al. ( 19991 )] or 



perpendicular |Patel etall (|l995 )]. depen ding on samp le processing details such as 
thermal history and shear starting time [Larso 3 (|l999h ]. At higher but still inter- 
mediate frequencies (where u < u> c , with u c the characteristic frequency of polymer 
chain relaxation dynamics), the preferred orientation is perpend icular for both PEP - 
PEE and PS-PI copolymers under high enough shear strain Koppi et al. ( 19921 ): 
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Figure 1: Three lamellar orientations (Parallel, Perpendicular, and Transverse) under 
shear flow. 



Patel et all f|l995h ; iMaring and Wiesnerl (|l997h ; iLeist et all f)l999h ] . At high frequen- 
cies (uj > u c ) the orientation selected is different for PEP-PEE (perpendicular) than 
PS-PI (parallel). 

No basic understanding exists about the mechanisms underlying the above com- 
plex phenomenology of orientation selection despite intense theoretical scrutiny in 
recent years. For the frequency range u < uj c that we are interested in here, the de- 
tailed relaxation dynamics of the polymer chains within each block are not expected 
to be important; thus one adopts a coarse-graine d, reduced desc ri ption in terms of the 
local monomer density as the order pa rameter Leibler ( 1980| ) ; lOhta and Kawasaki 
(198£jj); Fredrickson and Helfand ( 1987 )]. Most of the early analyses of shear align- 
ment of copolymers relied on thermal fluctuation effects near the transition poin t 
Jodt, and focused on the role of steady shears. A study by Cates and Milner ( 1989j ) 
indicates that in the vicinity of the order- disorder transition, steady shear would 
suppress critical fluctuations in an anisotropic manner, increase the tra nsition tem- 
peratu re, and favor the perpendicular orientation. Further extension bv iFredrickson 
( 1994 ) by incorporating viscosity contrast between the microphases has shown that 
the perpendicular alignment becomes prevalent for high shear rates, while the parallel 
one would be favored at low shear rates. 

Later consideration was given to situations that are not fluctuation dominated 
(such as well-aligned lamellar phases or defect structures), regarding both stabil- 
ity and defect dynamics. Stability differences between uniform parallel and perpen- 
dicula r structures subjected to steady shears have been found [Goulian and Milner 
( 1995f )] through the consideration of anisotropic viscosities in a uniaxial fluid. Re- 
cent m olecular dynamics studie s as well as h ydrodynamic analyses of the smectic A 
phase poddemann et all (J2004); [Qui (|2006j )] have shown an undulation instability 
of parallel lamellae, as well as a transition from fully ordered parallel to perpendic- 
ular phases for large enough shear ra te. Regardin g the e f fect of oscillatory shears 



ular pnases tor large enougn snear ra te. K egarding tne e l tect ot oscillatory shears , 
an analysis of secondary instabilities Drolet et al. ( 1999h : Chen and Vinalsl ( 20021 )] 



has shown that the extent of the stability region for the perpendicular orientation is 
always larger than that of the parallel direction, and as expected, both much larger 
than the transverse region. Imp ortantly, the role of viscosity difference between the 
polymer blocks (as introduced bv lFredricksonl ([19941 ) to address shear effects near the 
transition point) is found to be negligible for the stability of well-aligned lamellar 
structures, due to its weak coupling to long wavelength perturbations. In order to 
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address experimental phenomena related to dynamics of structure evolution and ori- 
entation selection, more recent theoretical efforts focused on the dynamic competition 
between coexisting phases of different orientations. Examples include the study of 
a grain boundary separating parallel and transverse lamellar domains under oscilla- 
tory shears that described the dependence of the grain boundary v elocity on shearing 
param eters such as frequency and amplitude |Huang et al. iHuang and Vinald 

( 200l ]. 

However, we are still far from accounting for existing experimental phenomenology 
on orientation selection, possibly because current approaches and models for block 
copolymer dynamics might not be adequate. It is important to note that block copoly- 
mer viscous response is not Newtonian even in the limit of vanishing frequency uo — > 0, 
since Newtonian response would result in the degeneracy of parallel an d perpendic- 
ular o rientations, contrary to experimental findings. In the theory of iFredrickson 
( 19941 ). the Newtonian assumption is used for individual microphases, with differ- 
ent Newtonian viscosities chosen for different monomers (blocks). We adopt here an 
alternative approach appropriate for flows on a scale much larger than the lamel- 
lar spacing and address the resulting deviation from Newtonian response in the low 
frequency limit. We introduce a constitutive law for the viscous stress tensor that ex- 
plicitly incorporates the uniaxial character of lam e llar phases in analogy with similar 
treatme nts of anisotropic flu ids Ericksenl (196Cf): Leslie! ( 196fih| and nematic liquid 
crystals jForster et all (|l 97lh: Martin et all (jl972f):lde Germes and Prostl iffooah]. As 
will be shown below, the resulting effective viscosity depends on lamellar orientation, 
in qualitative agreement with experiment results. 

The focus of our analysis is the competition among coexisting but differently ori- 
ented lamellar domains under oscillatory shears in a polycrystalline sample. We are 
primarily concerned with a simplified configuration involving two regions of parallel 
and perpendicular orientations, and use our results to address the selection mechanism 
in multi domain configurations. It will be argued below that domain competition, 
and the orientation selection mechanism it provides, could be determined by viscosity 
contrast between lamellar phases due to their different orientations with respect to 
the imposed shear. A hydrodynamic instability is predicted at the interface separat- 
ing parallel and perpendicular lamellae for certain ranges of material and shearing 
parameters. The resulting nonuniform secondary flows are argued to lead to the per- 
pendicular orientation in some ranges of parameters. Comparison of our results to 
existing experimental findings, as well as possible examination of our predictions, are 
also discussed. 



II. MODEL 

A. Governing equations 

In the range of low shear frequencies compared to the inverse of the polymer chain 
relaxation time, the phase behavior of a block copolymer is determined by an order 
parameter field ip representing the variation of local monomer density, and a velocity 
field v. Evolution of if) is governed by a time- dependent Ginzburg-Landau equation 

dip/dt + v • = -ASF /Sip, (1) 
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with T the coarse-grained free energy given by Leiblerl (|l980l ) and Ohta and Kawasaki 
( 1986f ). and A an Onsager kinetic coefficient. Equation f^Q) is coupled to the following 
equation governing local velocity field v = (v x ,v y ,v z ) for a block copolymer system 
under oscillatory shear flow: 



Re (dv/dt + v ■ Vv) = -Vp + V ■ a 



D 



(2) 



with the incompressibility condition V ■ v 
the Reynolds number defined as 



0. Here p is the pressure field, and Re is 



Re = pud 2 jf] 



(3) 



with u the shear frequency, d the thickness of copolymer system confined between 
shear planes, p the copolymer density, and rj a Newtonian viscosity. The coupling 
in Eqs. (0) and (J2J) is complex, especially when the stress tensor cr D in Eq. (J2J) 
depends on the concentration field ip, and when fluid inertia cannot be ignored (i.e., 
at nonzero Reynolds number). In order to simplify our analysis, we will assume that 
order parameter diffusion in Eq. (0) is negligible, so that ip is advected by the flow 
v. At the end of the analysis, we will discuss possible implications of the flow fields 
obtained on order parameter diffusion as given by Eq. ((TJ. 

Equation (J2J) has been made dimensionless by introducing a length scale d, a time 
scale to" 1 , and by rescaling pressure p — > pj (rju). In addition to a Newtonian viscous 
term, Eq. (J2J) includes a dissipative stress tensor cr D appropriate for a phase with 
uniaxial sy mmetry. We assume the consti tutive equation (in dimensionless form) due 
to Ericksen Ericksen (|l96nMLeslie1 as originally derived for anisotropic fluids 

and later applied to a nematic liquid crystal 



&ij = D i:j + aihihjh h hiD k i + a b& (hih k D jk + fijh k D 



ik i 



(4) 



where Dij = diVj + djVi = x,y, z), but with n 



a unit vector defining 



the slowly varying local normal of the lamellar phase. Eq. (J3J) is expected to be generic 
for any uniaxial phase, and thus should apply to the lamellar phases studied here by 
reason of symmetry. There are two independent viscosities ot-\ a nd (the Leslie's 
coefficients in the notation used by de Gennes and Prost (1993)), which have been 
rescaled as 

ax/rj and a 56 -> a 56 /r]. 



They o bey the relationship (due to the positivity of entropy production |Forster et al 
( 197l[ ): lMartin et all (|l972h ]) 



«56 > — 1 an d 0(x + 2a 56 > —1. 



(5) 



For a uniform lamellar phase of parallel orientation (n = (0, 0, 1) for the imposed 
shear flow along the y direction considered, see Fig. |TJ), the Navier Stokes equation 
can be written explicitly as 



Re (d t v x + VjdjV x ) 
Re (d t v y + VjdjVy) 
Re (d t v z + VjdjV z ) 



-d x p + d)v x + a 56 (d 2 z v x + d x d z v z ), 
-d y p + d 2 v y + a 56 (d 2 z v y + d y d z v z ), 
-d z p + (1 + a 56 )d 2 v z + 2(a t + a 56 )d 2 v z 



(6) 
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after substitution of Eq. (j^J) into Eq. (J2J). For a domain of perpendicular orientation 
(n = (1, 0, 0)) we have instead 



Re (d t v x + VjdjV x ) 
Re (d t v y + VjdjVy) 
Re (d t v z + VjdjV z ) 



-d x p + (1 + a 56 )djv x + 2(«i + a m )dlv x , 
-d y p + d]v y + a 56 (dlv y + d x d y v x ), 
-d z p + d)v z + a 56 (dlv z + d x d z v x ). 



(7) 



The Newtonian limit is recovered by setting ct\ = = 0. C ontrary to previous 
work o n block copo l ymers in the creeping fl o w ap proximation Goulian and Milner 
(1995J); Fredrickson ( 1994 ): Chen and Vinalsl (|2002l )]. we retain inertial terms in the 
above equations. Although small, they are significant in determining the instabilities 
discussed below. In the regime of small Re, both uniform parallel and perpendicular 
configurations are hydrodynamically stable as derived from Eqs. (JBJ) and ((Zj) with a 
procedure similar to that shown in Sec. IIII.I for the grain boundary configuration. 

Our assumption, Eq. (fl|). is consistent with experimental determinations of the 
loss modulus G" and dynamic viscosity 7/ (= G"/u) for uniform phases of different 
lamellar orientations. Assuming a shear flow along y (Fig. Q), we have in the creeping 
flow limit Re —>■ (in dimensional form) 

- Perpendicular (n = (1, 0, 0)): rj' = 77; 

- Parallel (n = (0, 0, 1)): 7/ = 77 + a 56 ; 

- Transverse (n = (0, 1,0)): 77' = 77 + a 56 + ^a tti (a = 7sin(u;£), with 7 the 

shear strain amplitude); when averaged over period T, it is (t/)t = V + «56 + 

9 



(1+7 2 ) 



i72ttl- 



When ai > and a 56 > 0, we obtain > r^ lle] > ^ erpendicular , in agreement 

with low u; experimental results of Koppi et al\ 1 19921 ) for PEP-PEE dib lock c opoly- 
mers, and also with the result of molecular dynamics simulations by Guol (2006) 
showing smaller viscosity of perpendicular lamellae compared to that of the parallel 
phase. On the other hand, w hen < 0, we have //[,., n .| < r?p„ rF „ ndi ■ , which is 
the case in PS-PI copolymers Chen and Kornfieldl ( 19981: iGupta et a/1 (fl995h ] . The 
viscosity coefficients cii and a<~,§ could be in principle measured in block copolymers 
under shear flows, with experimental setups p ossibly similar to those used in nematic 
liquid crystals [see, e.g., methods reviewed by de Gennes and Prost (199c 



B. Parallel/perpendicular configuration and viscosity con- 
trast 

We focus on coexisting lamellar domains of parallel and perpendicular orientations 
under oscillatory shear flows. In real samples these two types of domains may be 
separated by topological defects such as grain boundaries, dislocations, or disclina- 
tions. We study here the simplified configuration shown in Fig. comprising two 
fully ordered, three-dimensional lamellar domains that are identical except for their 
orientation and thickness. The perpendicular domain A is of thickness d,A, and the 
parallel domain B of ds = 1 — d& (in dimensionless form), both confined between a 
pair of shearing planes. The system is uniformly sheared along the y direction; the 
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Figure 2: A parallel/perpendicular configuration subjected to oscillatory shear flow. 



velocity field v is zero on the lower boundary plane z = and equal to v = 7 cos t y 
(or 70; d cos (u>t) y in dimensional form) on the upper boundary z — 1. 

Under the shear considered here, bulk configurations of either orientation (parallel 
or perpendicular) are linearly stable. However, we will show that the configuration 
of Fig. El can become unstable. Note that the effective viscosities of the two domains 
are different, as discussed in Sec. IA~1 Therefore this configuration is analogous to the 
case of two superposed Newtonian fluids o f different v iscosity which is known to be 
unstable under ste ady (plane Couette flow Yih ( 1967 ): Hooper ( 1985^ ]) or oscillatory 
|King et~al. (jl999f )] shears. In the present case, however, the viscosity contrast follows 
from the orientation dependence of our constitutive law (£Q). This contrast between 
the lamellar phases can cause interface instability, but with a dependence on shearing 
conditions and system parameters much more complicated than that of the Newtonian 
limit. 



III. HYDRODYNAMIC STABILITY ANALYSIS 

A. Base flow 

The base state for the configuration of Fig. |2]is a planar interface located at z = d A , 
separating two stable perpendicular (A) and parallel (B) regions. Under uniform shear 
the velocity fields are along the y direction: v Aj b = (0, V AtB , 0). In dimensionless form, 
the velocities are given by 

V A = 1Ze{a A smh[(l + i)(3 A z] je a ] , 

V B = TZe{a A (sinh[(l + i)f3 A d A ] cosh[(l + i)p B (z - d A )\ 

cosh[(l + i)(3 A d A ] sinh[(l + i){3 B (z - d A )) ) 7 e**} , (8) 

with 

a~ A x = sinh[(l + i)j3 A d A ] cosh[(l + i)/3s<i B ] + ^/rncosh.[(l + i)/3 A d A ] sinh[(l + i)/5s<i B ] 
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and the viscosity ratio m = ha/ Vb (with ha = 1 for the perpendicular region and 
fiB — 1 + «56 for a parallel domain). Here (3a and are the inverse Stokes layer 
thicknesses 

/ D„ \ 1/2 

Also for the pressure field, = Pb = Po- Note that these base state solutions are the 
same as those of two superposed Newtonian fluids w ith different viscosities \ia and 
jj,B under oscillatory Couette flow King et al. ( 1999f )]. 



B. Perturbation analysis 

For Newtonian fluids or in some viscoelastic models (e.g., Oldroyd-B or Maxwell 
fluids), the three-dimensional stability problem can be reduced to an effective two- 
dimensional one, with fluid stability under shear flows governed by Orr-Sommerfeld 
type equations for a single stream function describing two-dimensional disturbances. 
However, this is not the case discussed here and governed by Eqs. (J2J) and (J!J), as will 
be seen below. 

We expand both velocity and pressure fields into the base state given above and 
perturbations, 

vf' B = V A , B 5iy + uf' B (i = x, y, z), Pa,b =Po+ Pa,b- ( 10 ) 

These flow perturbations are accompanied by an undulation of interface, denoted as 
((x,y,t) (see Fig. |2j). The boundary conditions at z = cLa + ((%,y,t) include: 
Continuity of velocity 

v A = v*, (11) 

continuity of tangential stress 

{ [1 " (d x () 2 }ag - d x ( d y C C - d y C ag + (ag - a° x )d x (} B A = 0, (12) 
{ [1 " idyCYK ~ d *t d v( <& ~ d ^ < + (?Z ~ °yy)dyC} B A = 0, (13) 

and balance of normal stress 

{-P + <& ~ d x ( ag - dyC a° z } B A = -T'(dlC + d 2 y (), (14) 
where {}^ = {}b — {}a, &ij is the dissipative stress tensor, and 

r' = T/( V cod), (15) 
with T the interfacial tension. Also, the kinematic condition at the interface yields 

(<9 t + v B - V)C = vf. (16) 
Assume expansions of the form 

u t' B = ut ,B {(lx, Qy, z, t) exp [i(q x x + q y y)} , (17) 

Qx,q y 

where q x and q y are wave numbers in the x and y directions. Substituting (|T7|) into 
Eqs. © and (J7J), retaining terms up to first order in the perturbation amplitudes, 
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and eliminating the pressure, we obtain the following equations for perturbed velocity 
fields u^' B and u x ,B , which govern the system stability when combined with the above 
interfacial conditions: 

For the parallel region B (d,A < z < 1), 

Re [(d t + iq y V B ){dl - q 2 )u B - iq y (d 2 V B )uf] 

= (1 + a 56 )(d 2 - q 2 fu B - 2a iq 2 d 2 u B , (18) 



Re [d t (d 2 z - q 2 )u B + iq y (d 2 z - q 2 )(V B u B ) + 2q x q y (d z V B )u B ] 

= (d 2 z - q 2 ) 2 u B + a 56 (d 2 z - q 2 ) {d 2 z u B - iq x d z u B ) - 2iq x a 1 d i z u B , (19) 

(here q 2 = q 2 + q 2 ), while for the perpendicular region A (0 < z < (1a), 
Re [{dt + iq y V A ){dl - q 2 )uj - iq y (dlV A )uf] 

= (d 2 - q 2 fu A z - a 56 ql(d 2 - q 2 )uj + tq^ql - a 5e (d 2 - q 2 )}d z ut (20) 



Re [d t (d 2 - q 2 )u A + iq y (d 2 z - q 2 )(V A u A ) + 2q x q y {d z V A )uj] 

= (1 + a 56 )(d 2 z - q 2 ) 2 ui - 2a iq 2 M - q 2 y )ut (21) 

with rigid boundary conditions on the planes z = and z = 1 

tf(0) = ^(0) = 0.#(O) = 5,^(0) = 0, 

fif(l) = «f (1) = ^f(l) = d z u B (l) = 0. (22) 

The other velocity field Uy' B can be obtained from the incompressibility condition. 
Note that the form of Eqs. (fTH|) - (j21j) is similar to that of the Orr-Sommerfeld equa- 
tion for Newtonian fluids; however, in the above equations velocity fields are coupled 
(except for q x = 0), and thus the three-dimensional system here cannot be reduced 
to an effective two-dimensional one described by only one stream function. This 
irreducibility might be understood from the fact that parallel and perpendicular ori- 
entations can be distinguished only in three-dimensional space, and thus the results of 
the associated flows that are orientation dependent should be also three-dimensional. 
The solutions to the above problem can be found by writing 

uj' B (z,t)=e° t <t>f' B (z,t), 

u^ B (z,t)=e^^ B (z,t), (23) 
((t)=e° t h(t), 

with interfacial perturbation £ defined by ((x,y,t) = Ylq x q C{t) exp[i(q x x + QyV)]- 
When q y ^ 0, according to Floquet's theorem (p A ' B and h are periodic i n time t wit h 
period T (= 2n here) if the eigenvalue a (the Floquet exponent) is simple Yih f(l968h ] . 



since coefficients in Eqs. f!TH|) - (j2*T|) that are proportional to the base flow Va, b are 
periodic in t as shown in Eq. (jHJ). When q y = 0, 4>f' B and h are time independent, 
and a represents the perturbation growth rate. In either case, the real part of a 
determines the system stability. 
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C. Small Reynolds number limit 

The Newtonian viscosity rj is very large for block copolymers, resulting in small 
Reynolds numbers (Eq. (jSJ)). For a typical block copolymer of density p ~ 1 g cm -3 
and 7] ~ 10 4 — 10 6 P, Re/u ~ 10~ 4 — 10 -6 s for a thickness d ~ 1 cm. Within a 
reasonable range of frequencies uj of interest, we have Re <C 1. The stability problem 
can be now solved analytically by expanding around small Reynolds number, 

^; x B (z,t) = 0f;* o (;M) + Re ^ B x)1 (z,t) + 

h = h (t)+Reh 1 (t) + (24) 
a = a + Re oi + 

with cf) and h functions defined in Eq. (|2*3*|) . Note that the order of the rescaled 
interfacial tension T', defined in Eq. fTPoj) is to dependent, and needs to be addressed 
separately. In the limit of Re — > with V finite, we set T' = T = 0(1), while 
for small but finite values of Reynolds numbers, T' can be expressed as a power of 
Re. For a typical v alue of T ~ 1 dyne /cm (as estimated from the interfa cial energy 
between two blocks [Helfand and Tagami ( 1972h : Gido and Thomas! ( 1994h ]) and large 



7] (~ 10 4 — 10 6 P) appropriate for block copolymers, we have r'u; ~ 10~ 4 — 10~ 6 
s -1 . The intermediate range of u (around 1 s" 1 ) for typical experiments leads to 
r' = Re Fx = O(Re), whereas for much larger frequencies T'/Re <C 1. In the 
following, we present the solutions of zeroth order (containing To) and first order 
(containing Ti) in Re, which are accurate enough to determine the stability behavior 
for lamellar block copolymers with typical Re <C 1. 

C.l Re -> (with r = r = 0(1)) 

In this parameter range, we only need the zeroth-order solution for the velocity fields, 
given by (for q x ^ 0) 

0f o = B^e^ + BPe^ + B^e^ + B^e^, (25) 

<f>*o = C[ 0) e r ^ + Ci 0) e r2Z + Ci 0) e r3Z + C^e r ^ + ^d z ^ , (26) 

eft* = Af ] e aiZ + Af ] e a2Z + Af ] e a * z + A < £ ) e a4Z 

+ Df ] A\e SlZ + Df ] A' 2 e S2Z + D { 3 0) A' 3 e S3Z + A' 4 e S4Z , (27) 
# = Df ] e siz + DP e S2Z + e S3Z + Df V 4 *, (28) 



with 



bt = — \fi B + ui± [«i(2/i B + «!)] 1/2 | , (29) 
s 2 i = {nBq 2 + aiql ± ql [oti(2fi B + "i)] 1/2 } /fJ-B, (30) 

* = «* ( , 2 !; 4 , - 4] « (3D 



- q* - a^q* <& 
for i — 1, 2, 3, 4, and 

n = q, r 2 = -q, r 3 = qj ^/JT^, r 4 = -qj Ju~~b~, (32) 
ai = q, a 2 = -q, a 3 = \/ q 2 + a^q 2 , «4 = -a 3 . (33) 
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For q x = 0, we have = 0f o = 0, 0f o is also given by Eq. (|23)). and 

= (4 0) + 4 0) z)e 9Z + (-4 0) + 4°^)e" 92 - (34) 

By using the boundary conditions (J22j) . we can express the coefficients A34, -83°], 

Cg ], and D3 ] in terms of Afl, B±\, C[°l, and D^\. The remaining coefficients are 
obtained by solving a linear matrix equation 

X -AW = h®ho, (35) 

where A (0) = [Af , sj° >, cf, £>< 0) ] (j = 1,2) for g x ^ or [Af , ] for fe = 0, X 
isa8x8or4x4 constant matrix determined by interfacial conditions (JTTj) -(|14| ) . 
^0 — ho{t) denotes the Ot/i-order interface perturbation, and the matrix is a 
function of T and (d z V B ) , with ((9 2 Vb)o the gradient of Ot^i-order base flow V^: 

v b ] = [-Ma + S(z - d A )]^cost, (d z V B )o = £7 cost, (36) 

with 

5 = m\ , X = (d A + md B )~ l (m = fi A / fi B ). 
To lowest order, the kinematic equation (fT^j) yields 

d t h = -((To + iq y V B0 )h + <f>z {d A ), 

where V B o = Vg(z = d A ) = Ao^aT cost. From Eqs. (f2~K|) and (jSHJ), 4>fo(dA) can be 
expressed as 

4%>(dA) = [f£,i(qx, q y )r + f%(q x , q y )(d z V B ) ] h , (37) 

with fz 01 and f z02 complicated but known functions of wave numbers q x and q y (and 
also dependent on parameters a±, a 56 , and d A ), and hence 

d t h = [-cr + /Jo,i(?<d> 9y) r o] ^0 + [-W^Vbo + f^M^ 9»)( 9 *Vb)o] ^o- (38) 

The requirement of periodicity of /iq with time t gives the Ot/i-order Floquet exponent 
(note that terms proportional to V B0 and (d z V B ) are periodic in t) 

ff o = fzo,i(qx,q y )r . (39) 

Consequently, with initial value h (0), 

h {t) = h (0)exp{j dt[-tq y V m + f^ 2 (d z V B ) }} 

= h (0) exp [(-iq y d A + m/^ 2 )A 7sint] , (40) 

and the velocity fields are obtained from Eqs. (J2HJ), (J25j) - (J25|l . and all propor- 
tional to /io(£)- 
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C.2 Re < 1, r < 1, such that r'/ite = V x = 0(1) 

In this case, the results of first order expansion need to be addressed. When q x ^ 0, 
the solution to O(Re) yields 



= £ iB^e*" + B' iZ e biZ + B' i+4 z 2 e b ^ , (41) 
= E^ + ^0* 

i=l ^ 
4 4 4 

+ e c 'i^ nz + E + E ^' eM ' ( 42 ) 

i=l j=3 t=l 

4 4 

= E [^ 1} e a!2 + A a) ^i eS<2 + A'f'ze^ + E ^i+2^ 2 e aiZ 

i=l i=3 
4 

+ E We*** + D 'U* eSiZ + D 'Us^ SiZ ] (43) 
»=i 

4 

E WV«* + A'^ 2 + £>- +4 ^ 2 e Si2 + 4'e a < 2 ] , (44) 



»=i 

4 



i=l 



while for g x = 0, we get fa x = 4> xl = 0, 0f x given by Eq. (jHJ), and 

<jfa = (A? + z)e qz + + z)e- qz 

+ {A'l'z 2 + A%z 3 )e qz + {A'^'z 2 + A'i'z 3 )e- qz . (45) 



Here exponents 6j, r«, ai, Sj and coefficients A\ are the same as those in Eqs. 
(J33|), and coefficients S^, Sf, A' { , A", A'", D'^ and D" are complicated functions 
of 0th order solutions (A^, C^°\ D^). The unknown first-order coefficients 
(A« 5« C« L>W) can be determined by boundary conditions (|22jl and interfacial 
conditions (jll|) -(|14j ) . as in the above 0th order case. Similar to Eq. ()35|). the linear 
matrix equation governing lst-order coefficients (= cf \ £)j 1} ] for ^ 

0, or [Af\ Bf ] ] for q x = 0, with j = 1, 2) is 

X- = h^ho + h^h, (46) 

with matrices x and /i^ -* the same as those in Eq. (|35|). the first-order inter- 

face perturbation, and the matrix a complicated function of Ti, (<9 z Vb)o, and 
^(c^VbV Thus, the solution of Eq. (|4l)Jl takes the form 

A0-)=[ x - 1 h^]h o +[A^/h o ]h 1 , (47) 

with given by Eq. (jHSJ), and accordingly the first-order velocity fields can be 
calculated with the use of Eqs. pi j) -(j4l) j) . 

The O(Re) result of the kinematic condition (|T6*j) is given by 

dthi = -(o- + iq y V B o)hi - {<J\ + iq y V m )h + fa{d A ), (48) 
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where V B i = A 2 7sint is the first order base flow V B evaluated at interface z = d A , 
with 

A 2 = -(Xid A - \ d A /3), Ai = [(d 3 A + m 2 d 3 B )/3 + md A d B d]\l. 

The value of (p^dj^) is determined by solution (jUJ), and from Eq. (|4Tj) we find 

= ft,)ri + fl 2 (q x , q y )(d z V B ) + fl 3 (q x , q y )(d z V B ) 2 

+/*M(fc, ?»)^(^b)o] + [0fo(^)/M /ii, (49) 

with the Otft-order solution as in Eq. (J37|) . fzii — fzov anc ^ /^h (* = ^ 5 3, 4) 

obtained from first-order solution (|4T|) and functions of ai, a 56 , and e^. Substituting 
Eqs. ()49j). and (j2HJ) into (HHJ), and using the condition that hi is periodic in 

time, we obtain the first order Floquet exponent 

n = /5,i(ff., «&)ri + i*V/iU&, fc)> ( 5 °) 

and the corresponding interface perturbation 

= ^(0)^0(0) + /f li2 57sint + /f li3 5 2 7 2 (sin2t)/4 

-(ig y A 2 + 5/fi, 4 )7(l - cost)] h (t), (51) 

with ho(t) given in Eq. (It is convenient to choose h\(0) = at time t = so 

that the initial condition for ft is h(t = 0) = ho(0), independent of Re.) 

Equation ()50|) determines stability. The first term of the r.h.s. is proportional 
to rescaled surface tension and is always negative, indicating the stabilizing effect of 
surface tension. The second term, proportional to 7 2 , incorporates the effect of the 
imposed shear flow and tends to destabilize the planar boundary. Detailed results are 
shown in the next section. 

IV. RESULTS 

Given the results of the previous section, the stability of the boundary depends on 
the two viscosity coefficients a% and a^, on the domain thickness d A , the shear strain 
7, as well as the Reynolds number Re and the rescaled surface tension V (both uj 
dependent). Here we focus on three characteristic regimes of V (all with Re <C 1 
as appropriate for typical copolymers): V = 0(1), T'/Re = 0(1), and V J Re <C 1, 
corresponding to different ranges of shear frequencies. 

A. Re -> and V = 0(1) 

The first regime of interest is that of very small Reynolds number (Re — ► 0) and 
finite surface tension (r" = To = 0(1)), which corresponds to very low frequency uj 
according to the analysis at the beginning of Sec. IIII.llC~l The Floquet exponent (or 
the perturbation growth rate) is then well approximated by zeroth-order result a in 
Eq. (|39|). Our calculations give a < (i.e., f ^1(^x1%) < 0) for all wave numbers 
q x and q y \ thus, the system is always stable, resulting in the coexistence of parallel 
and perpendicular domains under shear flow. The stabilizing effect of surface tension 
dominates, and cr oc T as shown in Eq. (jSHJ). 
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d A /d B =l/2 



<k/d B =2 



Figure 3: Growth rate a/ Re as a function of wave numbers q x and q y , for at\ — 1, 
a 56 = —0.9, 7 = 1, Re — 5 x 1CT 4 , as well as (a) = 1/3 (with d^jds = 1/2), 
with maximum growth rate cr max = 3.4 x 10 _4 i?e found at q max = (0, ±0.89), and (b) 
dji = 2/3 (with dA/ds = 2), indicating that a < at all wavevectors. 

B. Re «l,r'« 1, and Y'/Re = 0(1) 

In an intermediate range of lo, V is of order O(Re) as discussed in Sec. 1111. lid The 
perturbation growth rate can be written as a = Re o"i with the first-order exponent <j\ 
given m Eq. flHJJ). Note that x = f% t < whereas the maximum value of ff x 3 can 
be positive depending on system parameters oti, a^, and dA, leading to a competition 
between the stabilizing effect of surface tension (proportional to T\ = V / Re) and the 
destabilizing influence of the imposed shear (proportional to 7 2 ). Note that Ti can 
also be expressed from Eqs. (fT5|) and (jH)) as 

r x = 1/We = [T/(pd 3 )]uj- 2 = 9cu- 2 , (52) 

with We the Weber number and 9 = T/(pd 3 ). Thus, a system would be more unstable 
for larger shear strain 7 and frequency to. 

To examine the onset of instability, we have carried out a numerical evaluation 
of o for a range of typical system parameters. For a typical copolymer system, 
p — 1 g cm -3 , T = 1 dyne/cm, 7] = 10 4 P, and d = 1 cm, with Reynolds number 
Re = (10 _4 s)a; and rescaled surface tension T' = (10 _4 s~ 1 )/u;. The two dimensionless 
viscosities are set as: cti = 1, and a 56 = —0.9 (corresponding to an effective viscosity 
Pb = 1/10 and thus a ratio m = PaI^b = 10) or a^Q = 9 (corresponding to ps = 10 
and m = 1/10). 

Figure 01 shows the growth rate cti = a /Re as a function of wave numbers (q x , q y ), 
for OJ56 = —0.9 and two different domain thicknesses = 1/3 (with dAjds = 1/2) 
and 2/3 (with c?a/<^b = 2). The most dangerous wave numbers are near q x = 0, as 
can be seen in Fig. Et- The figure shows that the maximum growth rate (cr max = 
3.4 x 10~ 4 i?e) occurs at g™ ax = and g™ ax = ±0.89. Results for the larger ratio 
dA/ds = 2 are shown in Fig. Eb, indicating that the system is stable. When 0:56 = 9, 
we obtain opposite stability results: a small ratio o?a/^b = 1/2 corresponds to a 
stable configuration, whereas instability occurs for dAjds = 2, with cr max also found 
at q x = 0. It is also interesting to note that at dAjd^ = 1, instability occurs for all 
values of (i.e., for all effective viscosity contrast). Figure 0] shows the maximum 
growth rate for 7 = 1 and 0.5, uj = 5 and 10 (s _1 ) (corresponding to Re = 5 x 10~ 4 and 
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1x10"* 




Figure 4: Growth rate a/Re versus q y at q x = for c/^/cJb = 1 and different values 
of 7, u;, and a 56 . 

1CT 3 ), and two different viscosities a 56 = —0.9 and 9. As the shear strain amplitude 
7 or frequency u decreases, the range of unstable wavenumbers also decreases. This 
observation motivates the analysis of long wave solutions presented in subsection D. 

The perturbed velocity fields are defined by Eq. The velocity associated 

with the most unstable wavevector (q x = 0, q y ^ 0) is given by 

uT x = and <%°* = <i> z e X p{<T 1BSX t), 

with (p z determined by Eqs. ([25]). (J53j) . (jUJ), and Results for (fi z at the most 

unstable wavenumbers are shown in Fig. |3] (at t — T) , and in Fig. 13 (at the interface 
2; = (Ia), for ai = 1, 7 = 1, i?e = 5 x 10~ 4 (with uj = 5 s _1 ), as well as for a variety 
of domain thickness ratios cIa/cIb and different viscosity coefficients = —0.9 and 
9. The corresponding interfacial perturbation h(t), defined in Eq. (|2Hjl . is shown in 
Fig. Although the velocity fields near the interface are sensitive to the viscosity 
contrast and thickness ratios (e.g, the temporal dependence of 4> z for fi B = 1/10 
(with a 56 = —0.9) and fiB = 10 (a 56 = 9) is out of phase at the interface, as shown 
in Fig. EJ), the qualitative results are the same: The instability develops around 
the interface, and relaxes into the perpendicular (A) and parallel (B) bulk regions. 
Importantly, the perturbation flow is directed along the z direction (and y, due to 
the incompressibility condition), and hence it is transverse to the parallel lamellae. 
If the monomer density order parameter were allowed to diffuse, this secondary flow 
would result in the distortion of parallel lamellae, but not perpendicular. 

C. i?e«l and V / Re <C 1 

In the limit V /Re 1, which corresponds to small surface tension V or large enough 
lu, the effect of the surface tension is negligible at least for solutions up to first order. 
This leads to T = Ti = in the solutions of Sec. Illl.l Thus cr = 0, and according to 
Eq. (|3U)l the first order growth rate is 

= ^V/iU^)- (53) 

The function /^ 3 (q , x,%) can be positive, and note that it is independent of shear 
parameters 7 and uj. The growth rate is proportional to 7 2 . 
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Figure 5: Spatial dependence of the velocity amplitude <p z at time t = T (= 2n), at 
the most unstable wave vectors (g™ ax = 0, g™ ax ) given by Figs. ElandlU All the curves 
here correspond to unstable configuration with o > 0, with parameters Re = 5 x 10~ 4 , 
«i = 1, and 7=1. (a) a 56 = —0.9 (hb = 1/10), for cIa/cIe = 1, 1/2 and 1/3; (b) 
«56 = 9 ([is = 10), for d^/ds = 1, 2, and 3. The locations of domain interface are 
indicated by dotted lines. 
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Figure 6: Temporal dependence of the velocity amplitude <p z over a period T at the 
interface z = d^, with dAjds = 1, and other parameters the same as those of Fig. 




Figure 7: Temporal dependence of the interface perturbation h over a period T, with 
parameters the same as those of Fig. El 
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Figure 8: Velocity amplitudes <p x and <p z (at time t = T) as a function of position 
z, at the most unstable wave numbers (g™ ax , <?™ ax ) = (3.2,6.95). The parameters are 
chosen as a 56 = —0.9, d A /d B = 1, 7 = 1, and i?e = 10~ 2 (corresponding to large 
u = 100 s- 1 ). 



A typical profile for the velocity functions <ft x and <p z is shown in Fig. |H| In 
contrast with the regime discussed in subsection B (T'/Re = Ti = Oil)), the most 
unstable wave numbers g™ ax and g™ ax here are both nonzero. As shown in Fig. |H] 
and unlike the limit discussed in subsection B, perturbed velocity fields develop along 
both x and z directions (i.e., 4> x ,(p z 7^ 0), leading to the distortion of both parallel 
and perpendicular regions. Therefore, although the parallel/perpendicular interface 
will move in response to the hydrodynamic instability and the resulting secondary 
flows, the direction of motion and hence the dominant lamellar orientation cannot be 
deduced from this analysis. 



D. Long wave solutions 

The calculations presented so far show that in the limit Re C 1, f < 1, with 
T'/Re = 0(1), the wave numbers associated with instability lie at q x = and small 
q y , as seen in Fig. |H1 Thus, we discuss next a long wave approximation to the stability 
analysis. We expand the solutions of Sec. 1111. C in powers of q (= q y here), and find 
that 

fli = fli = -foq 4 + 0(q 6 ), (54) 



with 

and 
where 

with 



fo = ^d 3 A d 3 B (l + a 56 d A ) > 0, (55) 



fl 3 = hq 2 + f2q i + 0(q e ), (56) 
h = ^^d 2 Ad 2 BOt 56 f u f 12 , (57) 



fn = (1 + a^d 2 A ) 2 + Aa m d A d B {l + a 56 d A ) > 
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Figure 9: Stability diagram of thickness ratio d A /d B versus viscosity contrast fi B 
(= m" 1 ). For the unstable regime with a > 0, the perpendicular phase is selected 
over the parallel one, while the stability leads to the coexistence of parallel and 
perpendicular orientations. 



J2 



(d A - d B ){d\ + d 2 B ) + d\a A m + 2d\\d A d 2 B + 2(d A - d B )]a 



m 



1 - d\f - 2d\(\ + d B f]a 



+ 2d 2 A [3(d A - d B ) 2 + 
+ 2d A [2(d A - d B f + 3d A d 2 B - 4(1 + d A )d B }a 56 , 



(58) 



and f 2 a complicated but known function of d A , a>i, and a 56 . Here A is defined as 



A = (d B — n B d A ) + An B d A d B > 0. 

Therefore, the first-order Floquet exponent ()50|) can be rewritten as 

1 



(59) 



i 



c2fB 2 

+ 2 d **i,37 



2 5 hi 1 



(9f u;- 2 - -5 2 f 2l 



2 V. 



(60) 



Thus, when fx > at small g we have o\ > for all 7 and a;. From the definition of fi, 
we note that stability is determined by the sign of 056/12 (Eqs. (jBTj) and (jSHJ)) which 
is itself a function of a$Q and d A only, and independent of shear parameters 7 and to. 
The calculated stability diagram of d A /d B vs. fi B (= l + ase) is shown in Fig. The 
diagram is symmetric with respect to d A /d B — > (d A /d B )~ l and fi B — > fj, B . (Note that 
at d A /d B = 1 (i.e., d A = 1/2) instability is found for all values of \i B , in agreement 
with the numerical results in subsection |BT| This diagram reveals the analog of the 
thin layer effect in the interfac ial insta b ility c aused by viscosity stratification of two 
superposed Newtonian fluids Hooper! ()l985| )]. When the thinner domain has the 
larger effective viscosity, instability occurs. Unlike the Newtonian case, here the 
effective viscosity contrast is caused by the orientation dependence of the dissipative 
part of the stress tensor. Note the for a polycrystalline sample, although the value 
of Q!56 (fi B ) would be determined by the specific copolymer considered, the thickness 
ratio d A /d B would vary from domain to domain. Thus, according to Fig. El the 
development of the instability would differ in different portions of a large sample. 

The largest perturbation growth rate can be obtained from Eq. (jBUjl . When 
Ofo^ 2 - \Pf2l 2 > 0, or 

20f 



,2, ,2 



7W < 



5 2 f 2 



(61) 
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we obtain 



Jl 4 2 

7 w , 



16(0/ o - 5 2 /27 2 ^ 2 /2) 

gmax = H^o-^M 2 ^) 1U - (62) 

These formulae show that both maximum growth rate and most unstable wavenumber 
increase with shear amplitude 7 and frequency u, in agreement with the numerical 
results of subsection B and Fig. 0] 



V. DISCUSSION 

The analysis given is purely of hydrodynamic nature and makes no reference to the 
response of the lamellar phases to the flows considered. The fully coupled problem is 
very complex, but the flow analysis conducted here can be used to argue indirectly 
about orientation selection. The flow perturbation u will advect the lamellae through 
the advection term u ■ V"0 in Eq. (fl}. with ip the monomer concentration. Parallel 
lamellae are marginal to velocity fields along the x and y directions since these flows 
are parallel to the planes of constant ip, but will be distorted by flows in the z direction 
(Fig. |2J). Conversely, perpendicular lamellae are unaffected by flows along either z 
or y, but distorted by those along the x direction. The instability mode given in 
Sec. IIV.B for V /Re = 0(1) which might be of most experimental relevance (if we 
estimate the order of surfac e tension T from the polym er-polymer interfacial energy 



estimate tne order 01 suriac e tension 1 irom tne polymer-polymer mtertaciai energy 
|Helfand and Taeam 3 (|l972h : iGido and Thomas! jSj], as discussed in Sec. 1111. C 



and hence is our focus here, is associated with secondary flows with u x = and u z 7^ 0; 
thus parallel lamellae are compressed or expanded, while the lamellar configuration in 
the perpendicular region remains unaffected due to the absence of modulation along 
its normal. Distortion of parallel lamellae would create a relative imbalance of free 
energy T in the two domains: .^Parallel > -^Perpendicular- This free energy imbalance 
would be relieved through the motion of the domain boundary towards the distorted 
parallel region. Therefore we would anticipate that a consequence of the shear flow 
and the resulting interfacial instabilities would be the growth of the perpendicular 
region at the expense of the parallel one. 

Therefore, the stability diagram of Fig. |U] can be used to indirectly address ori- 
entation selection, suggesting coexistence of parallel and perpendicular domains (in 
the hydrodynamically stable regime) or the selection of the perpendicular orientation 
(in the unstable regime). It would be interesting to examine an experimental system 
composed of only two domains of parallel and perpendicular orientations, for veri- 
fying our predictions such as the change of instability with domain thickness ratio 
cIa/ ds and viscosity contrast /xg, measuring the viscosity contrast from the location 
of the instability boundary, and further studying the domain evolution beyond the 
instability stage. 

Polycrystalline samples of the type present in all shear aligning experiments in- 
volve a distribution of grain sizes and orientations. We address next how the results 
just obtained may pr ovide a criterion fo r orientation selec t ion u nder certain circum- 
stances. Experiments |Gido et ~al. ( 1993i l: Qiao and WinevI (|2000h ] reveal the presence 
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of grain boundaries separating domains of different orientations, and we argue that 
the motion of these grain boundaries under the imposed shear affects the selection 
process. This mechanism is different than other suggestions in the literature involv- 
ing grain rotation, domain instabilities, or other effec ts of microscopic origin that are 
related to block architecture (looping and bridging) |Wu 

Earlier research has shown that boundaries of domains with a lamellar normal 
that has a component along the transverse direction will move, leading to a decrease 
in the size of the domain. The shear increases the free energy density of the transverse 
domain and originates diffusive monomer re distribution at the b oundary to reduce 
the ex tent of the phase of higher free energy Huang et al. ( 2003h ; Huang and Vinalsl 
( 20041 )]. Since the free energy of neither parallel nor perpendicular lamellae is af- 
fected by the shear ( at least in the low f requency range of u -C ay in which the 
Leibler Leibler ( 1980h ] or Ohta-Kawasaki Ohta and Kawasaki (1986)] free energies 
are a good approximation), one can generally expect grain boundaries to move to- 
ward the transverse phase. Our interest in this paper is therefore in possible physical 
mechanisms that would account for the motion of boundaries separating domains of 
parallel and perpendicular orientations. According to Fig. El if > (as might 
be appropriate, for example, for PEP-PEE diblocks) an initially large perpendicular 
domain (A) adjacent to a smaller parallel domain (B) (so that (1a/ ds is large) would 
grow even larger. Although our analysis does not hold beyond the linear stage of 
boundary deformation, it seems unlikely that any nonlinearity could saturate bound- 
ary distortion and lead to a stationary but corrugated boundary. Therefore we would 
predict that the perpendicular orientation will be selected for a 56 > 0. If, on the 
other hand, a 56 < (as would be appropriate, for example, for PS-PI diblocks), the 
situation is more complicated. Instability now occurs for o^/g^b small, leading to 
growth of the perpendicular domain and hence to an increase of the characteristic 
scale dA- To the extent that, in a sufficiently large system the boundary remains 
quasi-planar, the stability boundary in Fig. El would be reached. Once inside the sta- 
ble region, any remaining curved boundaries would be expected to relax to planarity 
(driven by excess free energy reduction), as the planar boundary would no longer be 
unstable under shear. Therefore, in the case of a 56 < 0, we would anticipate coexis- 
tence of parallel and perpendicular domains, or perhaps a dependence of the selected 
orientation on initial condition or sample history. 

We finally examine the dependence of the largest growth rate cr max on the shear 
amplitude 7 and angular frequency u near the onset of instability. Note that the 
stability boundaries of Fig. El are independent of both parameters, but near onset 
where a ~ O(Re) <C 1, experiments might detect an effective stability boundary 
located at the point on which l/o" max is of the order of the observation time of the 
experiment. From Eq. (j3J) we have 

a™ ax = a max /Re = a max r]/(pd 2 uj), 

and then given Eq. (|62|) we find that 



2 2/2 , 8,/2Wmax \ 160/ O 77<W , fiQ x 



with associated wavenumber 

„„... = ( 

Ofopd 



l^max \ 1/4 



o/ /4 . (64) 
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Figure 10: Lines of constant <T max (= 1CT 11 ) for (1a/ ds = 1. Symbols (circles and 
diamonds) are from numerical calculations, while the solid lines in the inset follow 
the scaling given by the long wave approximation in Eq. (J65j) . that is, 7a; 3 / 4 = const. 



Usually 7 2 a; ^> 8\f2\na mauX /(5 2 ffpd 2 ) for very small value of cx max , and thus we obtain 

(65) 



160/oTJ 



which scales as 7a; 3 / 4 for a given block copolymer. Therefore, the line of constant 
(^max is given by 7a; 3 / 4 = const. This result is consistent with our numerical evaluation 
of solutions in Sec. IIII.C 2. as shown in Fig. EH for two cases of = —0.9 (circles) 
and 9 (diamonds). Note that the data in the log-log plot of the inset is well fitted 
to a straight line with slope —3/4. For small enough cr max (e.g., = 10~ n as in Fig. 
llOp . which is the order of the inverse experimental observation time, this line can 
correspond to an effective stability boundary above which the instability becomes 
experimentally observable. 

To our knowledge, the only experimental determination in 7-cu space of the regions 
in which parallel or p erpendicular lamellae are sel e cted has been carr ied out for PS- 
PI block copolymers Marine and Wiesner (1997); Leist et al ( 1999f )]. It has been 
found that the line separating regions of perpendicular and parallel orientations was 
approximately given by ^uo = const. No results for copolymers with > (such 
as PEP-PEE) are available. Given that our predictions of effective boundary only 
apply to this case, it would be of interest to repeat the experiments for this type of 
copolymers. 

The above discussion corresponds to a range of shear frequency so that V /Re = 
(9(1). For very low frequencies so that Re — > while T' = 0(1), a parallel/perpendicular 
configuration is always stable due to the dominant effect of surface tension (as given 
in Sec. HV^). 

In combination with the stability of the lowest frequencies, our results in Fig. 
ITU1 indicate that the perpendicular orientation would be observed for large enough uj 
and 7. Otherwise, both parallel and perpendicular lamellae would coexist, and the 
selection between them would depend on experimental details such as quenched or 
annealed history of the sample and the starting time of the shear, a s found in PS-P I 
copolymer samples [Patel et al. (|l995h : iMaring and Wiesnerl (|l997t) : iLarsonl (|l999h ] 
and cannot be addressed by the stability analysis here. 
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VI. SUMMARY 



The assumption of a dissipative part of the stress tensor a? which is compatible with 
the uniaxial symmetry of a lamellar phase leads to an effective dynamic viscosity 
that depends on the orientation of the lamellae relative to the shear. We expect this 
functional form of of- to capture the low frequency and long wavelength response 
to the lamellar phase, without making reference to the microscopic origin of the vis- 
cosity coefficients. We have explored here the consequence of this assumption on a 
configuration comprising parallel and perpendicular domains. Experimental evidence 
suggests that these two orientations are prevalent in shear aligning experiments, and 
we believe that the type of rheology proposed here may contribute to our understand- 
ing of the orientation selected as a function of the parameters of the block and of the 
shear. 

In particular, we have shown that an oscillatory shear imposed on a block copoly- 
mer configuration comprising lamellar domains of parallel and perpendicular orien- 
tations can cause instability at the domain interface. The instability manifests itself 
by finite wavenumber undulations of the velocity field along the direction normal to 
parallel lamellae, which we argue would ultimately result in the growth of the per- 
pendicular region at the expense of the parallel one. Our results indicate that the 
instability, and the selection of the perpendicular orientation, occur at an interme- 
diate frequency range of small but finite value of Re, and depend on both viscosity 
contrast and domain thickness ratio. This instability is analogous to the thin layer 
effect in stratified fluids; that is, the system is unstable when the thinner domain is 
more viscous. On the other hand, at very low frequencies (Re — > 0), coexistence of 
parallel and perpendicular lamellae is found as implied by hydrodynamic stability. 
Also, in contrast to previous studies, the selection mechanism between parallel and 
perpendicular orientations introduced here is of dynamical nature, and an indirect 
consequence of the secondary flows generated by the hydrodynamic instability of the 
two-domain interface. It would be interesting to test our predictions experimentally 
in a test configuration of block copolymers as we have discussed here. Note that the 
frequency range studied here is uo < uo c in which polymer chains remain relaxed and 
hence the details of the individual blocks are not important. Thus our results should 
be independent of number and type of blocks in a copolymer. 
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FIGURE CAPTIONS 



FIG. 1. Three lamellar orientations (Parallel, Perpendicular, and Transverse) under 
shear flow. 

FIG. 2. A parallel/perpendicular configuration subjected to oscillatory shear flow. 

FIG. 3. Growth rate a / Re as a function of wave numbers q x and q y , for ai = 1, = 
—0.9, 7 = 1, Re = 5 x 10~ 4 , as well as (a) = 1/3 (with dA/ds = 1/2), with 
maximum growth rate cr max = 3.4 x 10- 4 i?e found at q max = (0,±0.89), and 
(b) d& = 2/3 (with d^/dg = 2), indicating that a < at all wavevectors. 

FIG. 4. Growth rate a /Re versus q y at q x = for g^a/g^b = 1 and different values of 7, 
co>, and a§Q. 

FIG. 5. Spatial dependence of the velocity amplitude <p z at time t = T (= 27r), at the 
most unstable wave vectors (g™ ax = 0, (?™ ax ) given by Figs. H3 and 0] All the 
curves here correspond to unstable configuration with a > 0, with parameters 
Re = 5 x 1CT 4 , ai = 1, and 7 = 1. (a) 0:56 = —0.9 (/ie = 1/10), for ^a/^b = 1, 
1/2 and 1/3; (b) a$Q = 9 (/xb = 10), for dA/ds = 1, 2, and 3. The locations of 
domain interface are indicated by dotted lines. 

FIG. 6. Temporal dependence of the velocity amplitude <ft z over a period T at the inter- 
face z = dAi with dA/ds = 1, and other parameters the same as those of Fig. 

El 

FIG. 7. Temporal dependence of the interface perturbation h over a period T, with 
parameters the same as those of Fig. |U1 

FIG. 8. Velocity amplitudes <p x and <p z (at time t = T) as a function of position z, at 
the most unstable wave numbers (g™ ax , g™ ax ) = (3.2,6.95). The parameters are 
chosen as a 56 = —0.9, dA/ds = 1, 7 = 1, and Re = 10 -2 (corresponding to 
large to = 100 s~ 4 ). 

FIG. 9. Stability diagram of thickness ratio c^a/^b versus viscosity contrast \xb (= m^ 1 ). 
For the unstable regime with o > 0, the perpendicular phase is selected over 
the parallel one, while the stability leads to the coexistence of parallel and 
perpendicular orientations. 

FIG. 10. Lines of constant cx max (= 10~ n ) for dA/ds = 1. Symbols (circles and diamonds) 
are from numerical calculations, while the solid lines in the inset follow the 
scaling given by the long wave approximation in Eq. ()65j) . that is, 7a; 3 / 4 = const. 
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